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ABSTRACT 

Employing the L'Huillier, Redish and Tandy (LRT) wave function formalism 
we develop a partially connected method for obtaining few body reductions of 
the many body problem in the LRT and Beneze, Redish and Sloan (BRS) formalisms. 
This metnod for systematically constructing fewer body models from the N-body 
LRT and BRS equations is termed the Dominant Partition Method (DPM) . The 
DPM maps the many body problem to a fewer body one using the criterion that 
the truncated formalism must be such that consistency with the full Schrodinger 
equation is preserved. The DPM is based on a class of new forms for the 
irreducible cluster potential, which is introduced in the LRT formalism. 
Connectiv ty is maintained with respect to all partitions containing a given 
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partition which Is referred to as the dominant partition. Degrees of freedom 
corresponding to the breakup of one or more of the clusters of the dominant 
partition are treated in a disconnected manner. This approach for simpli- 
fying the complicated BRS equations ir appropriate for physical problems 
where a few body reaction mechanism prevails. We also show that the dominant 
partition truncated form of the BRS equations may be obtained by distributing 
the residual interaction in the exit channel in a manner consistent with the 
dominant partition truncations of the irreducible cluster potential. 


THE DOMINANT PARTITION METHOD 


I. INTRODUCTION 

Connected Kernel Equations (CKE's) have enjoyed considerable prominence 

In the recent history of reaction theory. The equations due to Alt, C.rass- 
1 2 

berger and Sandhas (AGS), Bencze, Redlsh and Sloan (BRS) and Kourl, Levin 
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and Tobocman (KLT) are representative examples. In the formalisms of AGS, 

BRS, and KLT the many body scattering problem Is formulated In terms of a set 

of coupled Integral equations for the transition operators. These CKE's are 

4 

often viewed as extensions of the three-body formalism of Faddeev . 

Although the CKE's provide mathematically correct formulations of the 
N-body scattering problem, these equations have not Inspired extensive 
uasge In direct reaction analysis. The Distorted Wave Born Approximation 
continues to be the primary method employed In the analyses of direct reactions. 
Significantly, there exists In the community an understanding that many body 
effects should be Included In reaction analysis;^ however, the CKE's are 
not generally regarded as offering a viable approach for such Inclusions. 

Even In the three-body case, the Faddeev Equations are often regarded as use- 
ful for mathematical proofs but not as feasible for calculations^. 

The complicated nature of CKE's as well as an uncertainty about how 
the dynamics Is distributed In these equations have been important factors 
in limiting the role of CKE's In reaction analysis. Necessarily any use of 
CKE's must involve truncations. This compelling necessity for methods of 
truncating CKE's probably has contributed to their limited use. 

The application of the CKE's to nuclear and atomic systems is rendered 
difficult because of the number of equations involved and the fact that many 
channels are theoretically treated on an equal footing. We therefore consider 
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what simplification can be achieved by reducing the number of equations 
and/or channels^ In practice since certain channels may be ignored or 
treated phenomenologically it may not be necessary to preserve connected- 
ness in them. We call the formalism derived from relaxation of full 
connectedness in a set of CKE's a partially connected formalism . A par- 
tially connected approach has been advocated as a means of circumv^enting 
the difficulty Imposed by the coupling of ali rearrangement channels by 
Hahn and Watson in the three-body problem^. 

The choice of criteria (by which one truncates a CKE) Is an open 
question. A possible approach to simplifying these equations for some 
problems is to map the many body space into that of a fewer body problem. 

This approach will be useful in the case that the physics seems to be 
dominated by a few-body mechanism. One example is the deuteron-alpha 
scattering at energies below the threshold for breakup of the alpha. 

Notably, such a mapping does not destroy all of the many-body infor- 
mation which would be lost in the arbitrary imposition of a few-body model 
on a given many-body system. 

Actually, when one considers such a truncation it becomes clear that 
many of the CKE's do not lend themselves to such a reduction method. The 
AGS and KLT equations are notable examples. The explicit dependence on the 
number of particles as exhibited by the AGS equations and the dependence on 
the number of channels as exhibited b> the KLT equations tend to make the 
structures of these equations rather rigid. The AGS equations make explicit 
the number of particles through kernels which contain all subsystem transi- 
tion operators. The KLT equations are written for a fixed number of channels 
which structurally excludes the possibility of later dropping one of the 


channels. 
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In this article we develop a method of truncating the many-body BRS 

g 

equations to a fewer body problem. The method is developed in the LRT 
connected kernel wave function formalism, and Is similar in spirit to the 
Hahn-Watson reduction method^. This truncation is termed the Dominant 
Partition Method (DPM) . It maps the given many-body problem to a fewer 
body problem whose solutions satisfy the full Schrodinger equation. 

The equations obtained constitute a partially connected set, the discon- 
nectedness appearing in those channels which are not considered explicitly. 

In section II the LRT wave function formalism and the related irre- 
ducible cluster potential are reviewed. In section III the Dominant Par- 
tition Theorem is presented and in section IV this reduction method is 
applied to the BRS equations. In section IV it is also shown that the 
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reduced sel of BRS equations may be obtained via a distribution method . 

The summary and conclusion are presented in section V. 

II. THE L'HUILLIER, REDISH, TANDY WAVE FUNCTION FORMALISM 

In this section a set of coupled connected kernel equations for the 
wave function describing the scattering between many-body (Ni4) clusters 
is obtained. These equations are derived by using the BRS equations and 
the Green function for the system. The system under consideration has N 
distinguishable particles which interact via two-body potentials. A 
division of the N particles into n clusters, is termed an a^-partltion. 

The Greek alphabet is used to label two cluster partitions and the N- 
cluster partition is labelled 0. The partition Hamiltonians H^, residual 
interactions V and associated Green functions are defined by 


6 


H = H 4- V , 

(1) 

a o a 

V® = H - H , 

(2) 

a 

G^ = (Z - H 

(3) 


where is the sum of two body interactions internal to the a-partition 
and Z is the complex energy parameter (Z = E 4 ie) . The full and free 
Green functions are given by 


< (Z - H)"^ 

(4) 

= (Z - 

o 

(5) 


where H is the full N-particle Hamiltonian and is the total kinetic 
energy operator. 

Consider the N-body scattering problem initiated by incoming bound 

states of the two clusters comprising the partition 6. The full wave 
9 

function is 

fg = Jllm ieG$g, (6) 


where <I>g describes a relative motion plane wave times the internal bound 
state wave functions for the two clusters. The Green function is express- 
ible in terms of G„, 


Using that 

we obtain 
Using (7) we have 
Noting that 


G = G_ 4- GV^G-. 

P P 

<=6*8 " <**^’^*S 


''e- 


'V ~ 
3 


(1 4- GV^)<i>g. 

<= <*s- 




(7) 


( 8 ) 

(9) 

(10) 
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where is the transition operator^® we write (9) as 


"a - ^ 


( 11 ) 


Employing the BRS equation , 


= Vg + K^°G T 

3 ^ a o 


ao„ 


( 12 ) 


we obtain 




,oo^ 


(13) 


In (12) V- is the sum of two body interactions internal to 3 and external to 

p 

11 2 

a. The kernel K^° is the sum of all Weinberg graphs of connectivity o 

o 

which begin with any interaction and do not end with an interaction in a. 
Defining K°° ^ K^j using (3) and (10) we have 


= 'J’o + I <3 K G V . 
3 3 ^ o a o 3 

a 


(14) 


This integral equation has a completely connected kernel. The operator 
is the sum of all o-connected Weinberg graphs. 

Decomposition of the wave function into parts associated with the two 
cluster partitions of the N-body problssn is achieved by writing 


’b ■ I 


(y) 


(15) 


where 


; (Y) ^ 

3 


$-6 + G K G . 

3 /3 o Y o 3 


( 16 ) 
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V Vy 

The wave function only has outgoing waves of the y-type, that Is 

bound clusters of the ypartltlon or direct breakup from the Y~partltlon. 
Equation (16) may be written In a convenient form by considering again 
(9) and writing G In terms of G . We find 




(17) 


Using (9) and (10) we have 

Multiplying from the left by G^G ^ we have 

G g“^'1'. = G G~^4>- + G 
oyB o3B o 3 


(18) 


(19) 


Using (19) In (16) gives 


= (6 - - G K G g"^)<J>. + G K G g"^4'„. (20) 

3 yB oyoB B oyoyB 


We not Introduce the operator V which Is defined by 


KG ^ Y Q . 
Y o Y Y 


( 21 ) 


The operator is termed the Irreducible Y-connected potential . It is 

the sum of Y-connected graphs which become less than y-connected if the 

12 

rightmost interaction is removed . We have 


(y) _ 


,-l. 


= (6 - - G 7 G G-^)<I>„ + G-7 , 

yB oyyB B OyB 


( 22 ) 


and G 7 G G.^ on sb ill is the same as G 7o6. . This yields 
oyyB oBBy 


= 6 .(1 - G 7-)$^ + G 7 'J'_. 
3 yB o 3 B o y B 


( 23 ) 
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Noting that 

(1 - G^V^S “ 

which is most easily seen on examining the anti-cluster expansion for 
(section III) we obtain 

t/Y) = G 7 (25) 

3 o Y 6 

We write this in differential form as 

« - «o - 

which is reminiscent of the Faddeev three-body result. These are the LRT^ 
equations. Ihe relation of these equations to other N-body CKE’s is dis- 
cussed in detail in ref. 13. 

III. THE DOMINANT PARTITION TliEOREM 


The anti-cluster expansion for the kernel in the BRS equation has 


been previously given^^, 


m=2 (a=')a m m 

m 


N-1 


.3 


(27) 
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The N's ir the above equation are termed counting coefficients . They 

depend on both a and a . From (21) we have the anti-cluster expansion 

m 

for the irreducible cluster potential. 


N-1 


V = J I N(a,a )V G G 
o m a a a 

m=2 (o»)a m m 

m 


-1 


( 28 ) 


Summing the components of equation (26) we note the interesting 
result that 
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(E - H - I E )T = 0. 
a 


( 29 ) 


This property provides the underpinning for the truncation of the BRS 
equations that will be presented in this section. Insight into the method 
is afforded by the following theorem. 

Theorem I. For arbitrary N, if V is given exactly then 


I I'/b - VVg 
a 


Proof : 


where V is the full potential and 'F solves the Schrodinger 

p 

equation for the N-body system. 

From (28) we have 


^^0 = 1 1 I N(a,a )V a G~^'F.. 

^ n / \ ’maaaB 

a a m=2 (oJ>)a m m 

m 


(30) 


Interchanging sums and using 


G~^ = G~^ - V a 3 a 
a a a ’ m 

m 


(31) 


we have 


N-1 


I I I -V^ v/)» 

a m=2 a a(aa ) m m m 


m m 


(32) 


Noting that 


a a 

V =V - V =V - V , a a 
o a a m 

m 


(33) 




N-1 a 

II I N(a,a ) V (1-G V ”* + G V> . (34) 

ra=2 a 0 (aa ) “ ^m \ ®m ^ 

m m 


we obtain 
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We use the result , 


'1' = 5, + G V “ 

e ga 6 a 6 

m m 


to obtain 


's ^ 1 . <=a ’b- 


m=2 a o(=>a ) 
n m 


m m 


In the appendix we show that 


I N(a,a = 


o(=a ) 
m 


a 

„ m 
m 


where c = (-1) (m-1) ! Using this in (36) and again employing (35) we obtain 

m 

n-1 

I ^ '1'q = I c I V 'P . (38) 

„ a 6 m ^ a„ 6 

a m=/ a m 

m 


Now we employ the lemmas: 


Lemma 1: 


^ c, N-l 
c k 
k 


Lemma 2: 


N-l f V 

\ "n 4-1 = 1. N - 3, 
n=2 


(k) (k) 

The S„ , are Stirlirg numbers of the second kind. S„ - is the number of dis- 
N-1 N-l 

tinct ways of making k clusters out of (N-l) objects. From (38) we now 
obtain 


y y f = V 4'„. Q.E.D. (41) 

^ a 6 6 

This means that (29) with the a sum taken over the two cluster partitions 
is the Schrbdinger equation. This result motivates the consideration of 

truncations of ^ and/or V such that the corresponding summed version of (29) 

a 
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remains the full Schrddlnger equation. We present a method which satisfies 
this condition through the following three results. 

Remark ; For a given N, If V ^ Is given exactly the arbitrary sum (J 
over a subset of the two cluster partitions does not result In (E - - 

\ = 0 becoming the Schrddinger equation when the sum is over an arbitrary 

a 

set of partitons a. 

Proof : Consider the case N = 3 for which V = V . We label the possible 

a a 

o's as = (1)(23), = (2) (13) and = (3) (12). Note that 


j *'o % - J \ ’e- 

01.02 Oi.Oj 


^1"‘^2 


(42) 

(43) 


Remark ; For a given N and truncated arbitrarily then the sum over all 

T 

a does not result in (E-H - [ 7^=0 becoming the Schrodinger equation, 

^ ° all a 

where V is a truncated version of 7 . 

a o 


I ’roof ; Consider N - 4 and suppose we truncate by taking ~ ^ a' 


A OW 


y 'y. = y V 4<.. 

^ CT 8 _ a 3 


all a 


all a 


(44) 


r (2) 

From Lemma 1, ) V = S„ ;V and we have 
’ “ 0 N-1 

-0 ‘3 “ “3 


I = sf >V Q.E.D. 


all 0 


(45) 


We now introduce notation to represent a particular class of truncations 


of the irreducible cluster potential. The operator is defined as a trun- 
cation of V that terminates with a single a , 3 m N - 1. This operator 

0 m 

Imrolves the set of partitions a such that a, $ a This truncated irreducible 

K K Ix 

cluster potential is defined in terms of the anti-clus • expansion for (28) as 
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^ N(a,a )V G G (46) 

° n=2 (ao)a^^aJ " 

n m 

This class of truncated operators allows us to Introduce the Dominant Partition 

Theorem (DPT). The term dominant partition derives from the role played by a 

fixed partition a in the truncation of V and in limiting the sum on two cluster 
m a 

partitions. 

Theorem II . (Dominant Partition Theorem) 

For arbitrary N, and an arbitrary fixed partition a^ 

I F “ ¥ = V >1' 

a(=>a ) ^ 6 $ 

m 

where V solves the Schrbdinger equation for the ^^-body system (3<m<N-l) 

p 


Using the definition (46) gives 


I = I I I <=a ’'s- 

a(=»a ) o(^a ) n=2 (aa)a (^a ) 

m m n m 


We write (46) in a more convenient form by picking off the a term 

m 

a m -1 , 

o(^a ) a(=»a ) n =2 (a=>)a (aa ) ^n ^n ^ 

m m n m 

+ y N(o,a )V G g“^ 

- ^ X m a a a 8 
a (^a ) mm 

m 

Again employing (31) we obtain 

^ ^o“ *^8 “ I I N(c,a )V (1-G V^") 

a(>a ) ° ^ a^a ) n =2 ( 0 »)a (^a ) " ®n ®n 

m m n m 


+ N(a,a^)V^ (1-G V "‘)]Y.. 
m a a 8 

m m 
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Using (33) and (35) we have 

a ra-1 

I '»'« = V $ + I I I N(a,a^)V^ V V 

o(=*a ) ^ ^ ® cr(aa ) n=2 (a»)a fea ) ^n ®n 

m m n m 

(5 

^ i “a *6- 

a(»a ) mm 

tn 

Interchanging sums we have 

a m-1 

I 1' ” f„ - V„ ♦„ + y I i N(o,a )V G V 

a(ia ) “ e 6 8 a (=a ) " “n % ® 

m n m 

(5 

i »<"-V''a '^a 

a(=>a ) mm 

m 


Use of (37) and (35) yields 


1 = t y 1 c V + C V ] ¥_. 

^ - R ^ -^ na ma B 


o(;>a ) ® ^ n=2 a (sa ) 


We now use 


Lemma 3: 


y v^® = v^® 

, V a m-1 

a (-3a ) n 
n m 


A proof is given in the appendix. This yields 


m-1 , . / \ a 

y 7 “m ^ = [ y C (S^"^V + V ®) + C V ] f„, 

V 0 n m a m-1 m a ^ B 

a(=>a ) B n=2 m m 

m 


where we have used (33) (with o replaced by a ) to express V in (52) as 

n a 

a n 

V + V ®. 
a_ a„ 


Noting the results 


y C = 1 - C 

^ n m n 


( 55 ) 
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and 


m-1 

1 c 

n m-1 

Tl-Z 


(n) _ 


= 1 


(56) 


which follow from Lemma 2 we have 


a(oa ) 
m 


a 

7 ® VJ/ 

a 3 


= V V 


3 


Q.E.D. 


(57) 


Theorem II provides the basis for the DPM. It shows that we may truncate 

through the anti-cluster expansion by retaining only those partitions 

that contain a given dominant partition. Note that the full partition 

Green functions G are retained in (46). They are not projected on the 

a 

n 

Hilbert space corresponding to bound states of the dominant partition a^. 

The description of the breakup of these clusters is contained in these Green 
functions. 

The reduced problem is then solved by solutions to the original 
Schrbdlnger equation. This theorem provides a consistent means of reducing 
the many body problem in the LRT wave function formalism to a few body 
problem. 

Nuclear reactions are commonly analyzed in terms of a few body picture. 
For a given N there are ^ - 1 two cluster channels. Any realistic 

attempt to solve the many body problem cannot treat all of these channels on 
an equal footing. Moreover, it is reasonable to expect that in direct 
reactions the processes involved are not so extensive that all possible re- 
arrangement and inelastic processes must be Included. In many cases a 
realistic approach to many-body reaction theory will be afforded by sytem- 
atlcally building a few body models. 
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IV. THE DOMINANT PARTITION AND THE BRS EQUATION 

We now obtain dominant partition truncations of the BRS equations. 

This is accomplished by restricting the sum on two cluster partitions to 

the class defined by o(»a ), where a is taken to be the dominant partition. 

m m 

Correspondingly we Introduce the appropriate truncation of the BRS 

kernel by using the anti-cluster expansion (27) . Restricting the sum on 

two cluster partitions to those that contain a particular a^ terminates the 

anti-cluster expansion with that term explicitly involving a . We write 

m 

^m 

the truncated kernel as k„ , that is 


m 




m 

n=2 (oP)a^( 2 a^ 


(58) 


n 


With 6,a( a ) we write the truncated BRS equations as 
m 

= V® + T C t"“. 

“ o 0>am) ^ 


(59) 


I 


/ 


It will be recalled that in the derivation of the BRS equations that a 
crucial step was the democratic distribution of the residual interaction over 
all partitions. If we restrict the distribution to those partitions contain- 
ing a and proceed with the derivation, are the truncated equations obtained 
m 

the same as those we have termed the dominant partition truncated BRS 
equations? The answer is yes and provides the next theorem. 

Theorem III . The dominant partition truncated BRS equations (59) are 

obtained by distributing the residual interaction over 
the subset of all possible partitions containing a^ and pro- 
ceeding as in the derivation of the BRS equations (ref. 16). 

This theorem provides a satisfying degree of consistency in the reduction 
of the BRS equations. These considerations are displayed in Fig. 1. 
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Formal expression 


distribution 
over all 
partitions 


distribution 

over 

dominant partition 
set 


Partially Connected 
BRS-type Equations 


anti-cluster 

truncation 


BRS Equations: 
Completely 
connected 


anti-cluster 

truncation 


DP-BRS 

Equations 


Fig. 1 


Proof : From (39) and (53) 


i ''1 = 


d. J 

J 


I V/ = V^“. 

k) 


j ro 


Also note that 


Using (56) we obtain 


I V^, B 

r V d. m-1 ’ 


d . (=»a ) j 
J ® 


^ 3 • 

m 


''* ■ I I C V® , e = a . 


We use this in the definition of the transition operator 


= V® G g"^. 
+ a 


We obtain 
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m-1 


d.a 


■ J! I Cj ''d ‘’d 5! I C. V® G. T ■' , 8 


•=» a . 
m 
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Using the Llppmann identity ’ to transform the Born term gives the 

following equation for a new set of operators T which are equal to 
on the half-shell 


Or. o o d.a 

T = I I C 6 + ^ I C G T ^ , 6 a . 

j=2dj(-ajJ j=2d,(.-*aj ^ 


j 


This yields 


j6a ^ „g 


m-1 da 

+ y y C. Vj G - T J , a a , 3 
“ j-2djk^) i ^3 


m 


19 

We use the Yakubovskii cluster expansion to decompose the transition 
operator internal to partition d^ into pieces of different connectivities. 
We write 

m-1 


^15.- I . I 


i n=j (d ^)d n 
J n 


and 


T®° G = V® G . • J K®° G . 

■"“j ° ■‘j ■‘l n-1 (d,=)d„ ■>„ ° 

J n 


The anti-cluster truncation procedure is to replace Kj° G by < ™ for 

d o d 

n n 

N > m and to drop all other terms. 

We then get 

1 4 * OL 

fin. Q S 

- vf + I I C l I K ” T J . 
j=2 dj ^a^) J n=j (dj=>)d^ '‘n 


We Interchange the n and j sums realizing that m is the largest number of 
clusters that we can have in the limited space. We obtain 
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« o ®"1 an d -1 

T®“ - V« + I I J. I C V J G , 

“ n=2 d (oa ) ^ n j=2 d (?d ) ^ 

n m j n 


where (60) has been used. Employing 
Lemma 6:^^ ^ 


J-2 d Pd^) J 


yields 


m-1 a da 

T^“ = + y y 6„ t " 

a \ Bd 2n 

n=2 d (^a ) n 
n m 


and 


a 

Sot ,,S r ro -C^ot 

T = V*^ + ) T 

a V pQ 

o(pa ) 

IQ 


6, asa . Q.E.D. 
m 


V. SUMMARY AND CONCLUSIONS 

We have developed a generalization of the Hahn and Watson’s^ 'partially- 
connected" strategy appropriate for constructing r.-cluster models for N-body 
problems where it is intended that n << N. The cases treated are those in 
which the only channels treated explicitly are obtained by combining the clus- 
ters of an n-cluoter "dominant" partition, a^. We obtain n-body equations 
of the BRS type for transition operators and of the LRT type for wave func- 
tions. The resulting equations are connected in the degrees of freedom cor- 
responding to the relative motion of the clusters of a^, but not in those 
Internal to a single cluster of a^. Following Hahn and Watson we assume that 
these degrees of freedom are to be handled in some manner different from oper- 
ator Integral equations (e.g., by statistical or phenomenological methods). 

Our main result is that the BRS and LRT equations for the small number 
of clusters is in fact exact if the subsystem Green functions are put in 
from some other source. This means that no incoming waves associated with 
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channels breaking the cluster of a are to be admitted. Furthermore, It 

n 

does not matter whether one obtains the partially-connected equations by 

truncating the anti-cluster expansion for the kernel or by only distributing 

the residual potential over the appropriate limited set of partitions. The 

same equations are obtained by both procedures. 

A specific example where a procedure such as described here" may be relevant 

Is In the slx-nucleon problem where the Initial channel Is a low energy 

4 20 

(E < 20 MeV) deuteron Incident on a He nucleus. As Is well known this Is 

well described as a three-body problem. Here, our dominant partition would 

be a^ = (n) (p) (nnpp) where the effects of exchange are Ignored. The three- 

body equation would fall out Immediately upon approximating the Green function 

4 

G by its part having the He pole. This approximation would yield real 

^3 4 

effective nucleon- He interactions. 

More general results, including the appearance of more complete effective 

Interactions, can be obtained in a number of ways, the simplest of which is 

the introduction of projection operators at the Green function G . The 

4 ■ ^3 

part corresponding to everything but the He pole are then solved formally A la 
21 

Feshbach . This leads to the appearance of generalized optical potentials 

22 

as effective interactions plus the well-known effective three-body force. 
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APPENDIX 

We now prove results which were employed in section III on the DPT. These 


results are Eq. (31) and Lemma 3. According to (37) 

I \Ha,a )V° = C V ", 

a(Da ) ™ “ 

in 


(A.l) 


where = (-1) (m-1) ! . An operator such as the T matrix for a given partition 
can be written in terms of a cluster decomposition 


~ ^ I ^a 2b f^^b » 


(A. 2) 


m=n b n m 
m 


where [T] is the part of T with connectivity a and 6^ is one provieded 

n ri~ m 

a contains b otherwise it is zero. Inverting the above expansion we obtain 
n m 

the [T] in terms of the T . The inverted expansion is termed an anti- 
m “ 

cluster expansion. In this expansion the numerical coefficients, N (the counting 
coefficients) occur naturally. From the cluster expansion for T„ we have 


[T], = I I N(b a ) T„ 

^n m=n (b 3)a ^ ^ 

n m 


(A. 3) 


The existence of the N(a ,a ) follows from the fact that the matrix <5 ^ -jv, is 

11 III 
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invertible . Moreover using (A. 2) and (A. 3) we easily note that 

k 

m=j a(cc.) m-k ik 

We develop a recursion relation among the counting coefficients by following 
Ref. 24. We write (A. 4) as 


Z Z +«c.=a >«a =>b, = «c.b, 

m=j a^ j- m j“ in m~ k j k 

where we have used 6 +6 = 1. We obtain 

c.^a c.sa 

J- m J- m 

k 

I I N(c.,a )fi 5 u " <5 u 

J m c,^a a =>b, c,b, 
m- j a J- ra nr k j k 


(A. 5) 


(A. 6) 


24 

since N(c.,a ) = 0 for c.^a . 

J m J m 

Employing (A. 6) we easily obtain a useful relationship among the counting 

coefficients which allow?: us to determine N(c,,a ) in terms of the N(c.,a,), 

J m j 3 

N(c.,a..,) N(c. ,a From (A. 6) we have 

3 3+1 3 m-1 


k-1 

y N(c.,a, ) 6 6 , + T y N(c.,a ) 6^ ^ 

^ i K c,3a, a, 2 b, ^ . i m c 2 

j- k k m=j a^ 3 


which yields 


.a (A.7) 

“ m m~ k 


k-1 


N(c ,b^) -- 1 I me 

j m= -i J in Tu iC 


(A. 8) 


m=3 a, 


m 


We return to (A.l) and employ an inductive argument. We note that for 
m = 2,3 that 


I N(a,a_)/ - C- V 


c 

I N(a,a = C V 


Assuming that 


we show that 


I N(a.aj^)/ = V 




a 


Consider 


I - I I J "(‘’■•’J «b 

a o m=*2 b mm k+1 

m 


v". 


where we have used (A. 8). Using the assumption (A. 11) we have that 

I = - I I v”” 5^ 


m= 


=2 b 


j -'a, , . 

ur k+1 


(A. 9) 
(A. 10) 

(A. 11) 

(A. 12) 

(A. 13) 


(A. 14) 


m 


We write this as 


I - - I I 


"■2 


(A. 15) 


Note that 
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a (=>a ) 
n m 


in m-1 


(A. 16) 


, (n) 


where the S'**' are Stirling numbers of the second kind. This result follows 
m 

almost immediately from Lemma 3 which will be proven. Employing (A. 16) we 
find 

I N(c.a^j)v'’ - - I C^(S« - S« (A.17) 


Fiaally we write 

I = -CT C„ S« - 

a m=2 m=2 


and use the fact that 


16 


N-1 


I C = 1 

tn m N-1 

\+l 


m’ 


to obtain 


(A. 18) 


(A. 19) 


I N(o,a-^^)v" = V 

a 


We complete the discussion by proving Lemma 3; 


Q.E.D. 


I 


m _ „(n) 


m 


a (=>a ) 
n m 


n 


m-1 


Proofs Consider a particular pairwise interaction with j external to a^. 

This interaction will be internal to some number of the a clusters that include 

n 

a . We determine how many by considering an m-cluster system. In all of the 
m 

n clusters the particles of the pair j will have to be together so we construct 

n clusters out of the m clusters by initially joining the two particles of 

the pair j into a single 'particle'. Now the remaining (m-1) clusters may be 

joined into n clusters in any way. This number is This is Independent 

m— i 

of the pair index so every pair external to a^ will appear this number of 


times. This yields Lemma 3. 
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